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Abstract 

I A new exactly solvable model for the evolution of relativistic kinetic 

' system interacting with an internal stochastic reservoir under the influence 

of a gravitational background expansion is established. This model of self- 
interaction is based on the relativistic kinetic equation for the distribution 
^ ' function defined in the extended phase space. The supplementary degree 

\^ • of freedom is described by the scalar stochastic variable (Langevin source), 

which is considered to be the constructive element of the effective one- 
^ _ particle force. The expansion of the Universe is shown to be accelerated 

' for the suitable choice of the non-minimal self-interaction force. 

(N 

O ■ PACS numbers: 47.75.+f, 04.30.-w, 05.20.Dd, 05.70.Ln 

o 

O 1 Introduction 

o 

^1 Observations of supernovae at high redshift show that the expansion of the Uni- 

^1 verse is accelerating. To explain this fact in the framework of Friedmann - 

Lemaitre - Robertson - Walker (FLRW) model, one has to require that a cos- 
mic medium is characterized by negative pressure jU Ej. One of the possible 
explanations is the existence of dark energy (see, e.g., the reviews There 
5^ I exists a number of dark energy candidates, the best known are a cosmological 

constant and different quintessence scenarios ( El El 13 UHl UH 112] ) • Negative 
pressure may also be the consequence of self-interactions in gas models of the Uni- 
verse [ini CH Uni Cni CZ! • in particular, an "antifrictional" force, self-consistently 
exerted on the particles of the cosmic substratum, was shown to provide an al- 
ternative explanation for an accelerated expansion of the Universe |14[ IT^ . This 
approach relies on the fact that the cosmological principle is compatible with the 
existence of a certain class of (hypothetical) microscopic one-particle forces, which 
manifest themselves as "source" terms in the macroscopic perfect fluid balance 
equations. These sources can be mapped on an effective negative pressure of the 
cosmic medium. 

It is worth pointing that effective self-interaction forces can be regarded as a 
specific non-minimal coupling of the cosmic gas to the Ricci scalar, Ricci tensor 
and Riemann tensor. Generally, a force which explicitly depends on curvature 
quantities describes a coupling of matter with the space-time curvature which goes 
beyond Einstein's theory. However, mapping the non-minimal interaction on an 
imperfect fluid degree of freedom admits a self-consistent treatment on the basis 
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of general relativity. This may be considered as a gas dynamical counterpart 
to the non-minimal couplings of scalar fields or those of higher-order gravity 
theories. A (non-minimal) fluid interaction is designed so that it results in the 
desired cosmic evolution. Designing the coupling for description of a specific 
dynamics has already been used earlier for interacting two-component models 
[H]. Here this idea is applied to the case of a one-component fluid, which is 
self-consistently coupled with the Riemann tensor. As a characteristic feature 
of this approach, Hubble rate and deceleration parameter explicitly enter the 
microscopic djTiamics, giving rise to a self-consistent coupling of the latter to the 
gravitational field equations. 

The paper is organized as follows. In section 2 the formalism of one-particle 
distribution function in the extended phase space associated with a homogeneous 
and isotropic, spatially fiat Universe is established. In section 3 macroscopic 
properties of the ultrarelativistic kinetic system are discussed, an effective stress- 
energy tensor is obtained, and the conditions of accelerated expansion of the 
Universe for power-law and exponential scenarios are investigated in detail. A 
brief summary is given in section 4. 

2 Formalism of One-Particle Distribution 
Function in the Extended Phase Space 

2.1 Kinetic Equation 

The idea of the phase space extension, based on the covariant formalism of Car- 
tan's differentiation and integration, was proposed in ^SJ E]. Along a line of 
development of the relativistic kinetic theory (see, [201 1^ 122] and references 
therein) a numerous applications of that formalism have been elaborated by dif- 
ferent authors (see, e.g., |23l 1211 123 1201 12Z1 ) • 

The simplest kinetic system with supplementary degree of freedom can be 
described in terms of the 8-1-1 -dimensional scalar distribution function 



$ = $(x\/,u;), (1) 



depending on coordinates x*, particle momentum four- vector , and on the ran- 
dom scalar variable a;, which is called random Langevin's source. The distribution 
function satisfies the kinetic equation 

where T^i are the Christoffel symbols, associated with the background metric gik- 
Characteristics equations, corresponding to (2), form three subsystems: 

*: + ri,/^ = f. §^ = w. if=fL. (3) 

as as as as mc 

The first equation in (3) is the well-known equation of general relativistic particle 
dynamics under the influence of the force ,u^. The second equation is 



2 



the evolutionary equation for the scalar random variable cu, which is modeling 
the stochastic influence of the environment on the particle. Using (3), one can 
write the rate of evolution of the distribution function $ in the form 



ds 



dp^ duo 



(4) 



The moments of the distribution function (1) can be obtained by averaging over 
five-dimensional statistical ensemble: 



N\x) 



dPduO ■ $ • p\ T\x) = / dPdu ■ $ ■ j9>'= . 



(5) 



Denoting the term / rfo; ■ $ by f{x,p), one obtains from (5) the standard formulas 
for the particle number density vector N'^{x) and for the stress - energy tensor 
T^^{x). The standard definition for the entropy flux vector: 



S\x) = -ksc / dPduj - p 



In h^<i> 



(6) 



is used. The transport equations for the particle number, stress-energy and for 
the entropy have the form, respectively: 



ViN\x) = , VkT'''{x) = mc / dPdu ■ $ ■ F 



dp^ duo 



aix) = ViS' = kemc^ / dPduo ■ $ 



(7) 



(8) 



2.2 Cosmological Background 

The spatially fiat FLRW solution of Einstein's equations is considered: 



ds' = c'df - a\t) {dx^f + {dx' f + {dx^ f 



a{£ + W) 



8nG 



(9) 



(10) 



3 ' " 3 

Here S and V are the effective energy density and effective pressure, respectively. 
These quantities are known to be the eigenvalues of the conserved effective stress- 
energy tensor 



ik 

(eflf) 



(eff) 



0. 



The velocity four- vector is equal to = Sq, as usual. 



2.3 Force Field Modeling 

The motion of Brownian particle in the framework of classical dynamics can be 
modeled by the force three-vector 



j?= -A(v -[/) + <. 



:i2) 
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The first term of tfiis force is tlie Stokes friction force, it vanishes when the 
particle velocity three-vector v coincides with the medium flow three-vector U . 
The second term is the stochastic Langevin force with random three-vector ^. 
The covariant generalization of the Stokes and the Langevin forces is well-known: 

-^stokes = >^{^k- ■^i^) ^'^ ' ^Langevin = >^ {^k ' ^'^ ' ^^^^ 

The projector with respect to particle four-momentum in the parentheses provides 
the forces to be orthogonal to p*: F'^pi = 0. 

In this paper we consider the force four-vector 

F\x\p\u;) = ypi)Sl-p'p>] , (14) 

mc^ L J 

where q is considered to be a function of cosmological time and is a subject of 
modeling. It was shown in [21] that the existence of the force with such a structure 
is compatible with the symmetry requirement 

£(;gik = 2^ (^gik - "^^^j ' (15) 

formulated in terms of Lie derivative. The corresponding time-like vector is 
shown to exist for FLRW space-time, and the explicit example of the so-called 
generalized equilibrium states has been found in [221 ■ Starting from this fact we 
have studied in [HI [T7j the consequences of the appearance of the antifriction 
force. The force (14) gives the Stokes force when ujqm = X and the Langevin 
force when ujqmU^ = k^^. 

The function Ti is modeled as follows: 

H(x^p^a;) = ^^(//7,). (16) 

The coefficients q and x ^-^e considered to be scalar functions, depending on the 
Ricci scalar R, on the Hubble parameter H (which is proportional to the scalar 
of expansion H = ^VkU^) and on the scalar R = RikU'^U^ . The terms and Ti 
are linear in dimensionless random variable cu, i.e., one can indicate a stochastic 
Langevin's source uo as multiplicative one jSU]. The function Ti in the form (16) 
guarantees that the trivial value oj = is a singular solution of the dynamic 
equation (3). Since 



dF' m {p% 



9p* du mc^ 
the entropy production scalar a (8) does not vanish. 



+ — = ^^-^ix-^q^)^o, (17) 



2.4 Solution to the Characteristics Equations 

Supposing that in the FLRW space-time (9) all the macroscopic functions depend 
on time only, and using the consequence of the characteristics equations 

f = (18) 

as mc 
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one can rewrite the dynamic equations in terms of t 

Pa = -qujpa, io^x^- (19) 

The dot denotes the derivative with respect to time, and a = 1,2,3. The solutions 
of (19) are the following: 

Va = C«e-"-^(*'*°), J{t, to) = f dtq{t)I{t, to) , (20) 

Jto 

u; = n-I{t,to), /(t,to) = eX''*^(^ (21) 

By definition J{to,to) = and /(to,^o) = 1, providing the relations Pa{to) = Ca 
and cu{to) = ft. Using the solution for p^ and the normalization condition p^pi = 
m^c^, one can obtain the p^ — po component of the particle momentum: 

p°{t) = ^Jm^c^ + a-2(i)C2e-2f^^(*'*o) . (22) 
Here C'^ = + C2 + C^, and the distribution function $ takes the form 



$ = $o(C^ to)e'''"^'''^^I-\t, to)S{^{p,p) - m?c^) . (23) 
The transformation of the generalized volume element gives the formula 

[ dPdu .${•..}=/ , dCdCdCs ^^^^ 

J J a%t)Jm^c^ + a-2(^)C2e-2^''(*'*o) 



3 Macroscopic Properties of the Kinetic System 

3.1 The Structure of Macroscopic Moments 

Particle number is a conserved quantity, thus, one obtains 



N'^a^t) = const, N'{t) = N{to) 



q(^o) "'^ 
a{t) 



Sh ■ (25) 



The stress-energy tensor takes the form 
nit)- I , "'"'"j!'^, $o(C^Q,to)•r^(C^J(Mo)), (26) 



where 



Tl{C\J{t,to))^a-^C^e-^'''diag\[mca\C\-'e''']' + l, ~, . (27) 

In general case the stress-energy tensor can not be represented in the analytic 
form, thus, for the sake of simplicity the ultrarelativistic limit is used to illus- 
trate the idea. If the term, containing m?(? in the function (27), is negligible in 
comparison with the second term, and if the initial distribution function $o is 
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multiplicative, i.e., $o(C^,l^,^o) = /o(C^,to) " ^(^i^o)) one can obtain the stress 
- energy tensor in the form: 



mt) = A{t,to)w{to) 

Here 







a{t) 
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A{t,to) = J dn^o{^,to)e-^'^^'''°^ = {e-^-^) (29) 
is a statistical factor. The energy density scalar 



W{to)^J^^\C\fo{C',to) (30) 

is known to be equal to 

W^bosons(^o) = YS/JS^"^ ^^^^ ' ^fermions(^o) = 3Q/^3g3 -^ ^^^^ ^"^^"^ 

for ultrarclativistic (massless) bosons and for ultrarelativistic (massless) fermions. 
The entropy production scalar 

a(t) = kBcN(t) J dmoi^, to)(x - ^Q^I) (32) 
happens to be proportional to particle density scalar N{t). 

3.2 Properties of the Averaged Macroscopic Moments 

Let us suppose that the distribution over ft can be described by the Gaussian 
function 

^j,(n^to) = ^-e-^ , (33) 

where D is a dispersion parameter of this distribution at the moment to. For this 
case the direct calculations give 

a{t) = kBCxN{t) (34) 

and 

A(i,to) = e^^''^'(*'*°). (35) 

Thus, the properties of the macroscopic moments of the distribution function are 
predetermined by the properties of the functions J{t,to) (20) and I{t,to) (21). 

3.3 Effective Stress-Energy Tensor 

The relation (7) demonstrates that the energy and the momentum of the kinetic 
system do not conserve, since the balance equation has the source term in the 
right-hand side. Calculating this source term using the expression (14) for the 
force one obtains 

V.T^'^ = gl{t, to) / -7-T===T^r^[-'^'^^ - P'^U'^P'^^ ■ (36) 
ca%t)^'m?c^ + a-^{t)C^e-^^J 
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For the ultrarelativistic model with multiphcative distribution function this ex- 
pression is simphfied: 



a{t) 



-n.j 



|C|/o(C2,to)- (37) 



For the Gaussian function (33) the integrals in the balance equation take the form 



VkT"'^SlqI{t,to)W{to) 



a{t) 



(38) 



For i = 1,2,3 these equations are trivial, the only informative equation is the one 
for i = 0. This equation (the so-called scalar balance equation) may be written 
as follows: 



a{to 



a{t) 



J{tM)e-'^''''\ (39) 



W{t) + ?>H[W{t) + P{t)\ = ql{t, to)W{to) 

where W{t) — T°°. The equation (39) takes the form of conservation law 

W + m{W + 7^) = , (40) 
if we introduce an effective pressure V = P + 11, where 

a(to) 



n 



-qI{t,to)W{to) 



a{t) 



6H{t) 



(41) 



is considered to be non- Pascal pressure. The reconstruction of the effective stress- 
energy tensor (11) is possible if we put 



E = W{t) = W{to) 



a{to) 
a{t) 



1 4 



e 4 



V = -Wito) 



a{to) 


4 2 

e 4 


[ «(^) J 





l-qI{t,to) 



2H{t) 



J{t,to) 



(42) 



(43) 



The formulas (42) and (43) manifest the following feature. The coefficient 

describes the decreasing of the average energy and pressure, which corresponds 
to the standard law of evolution of ultrarelatistic FLRW model. The coefficient 

-D t2 

e~'' describes the increase of the mentioned quantities due to the gas (fluid) 
stochastic self-interaction. The interplay between these two processes could clar- 
ify the question whether £ is the increasing or decreasing quantity. The sign 
of V can be both positive or negative during the different time intervals. The 
expansion of the Universe happens to be accelerated (a > 0)(see, (10)), when 



£: + 37^ = 2vr + 3n < . 



It is possible when 



q{t)I{t,to) 



4H{t) 

Let us discuss this inequality in detail. 



J{t,to) > 1. 



(44) 

(45) 
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3.4 Modeling of the q and x Functions 

The functions q and x can be modeled using R, R and H functions: 



i? = -6 



d 

- + 




)] 


a 










-3 


a 







a 

a 



(46) 
(47) 



Taking into account tliese relations, q and x can be rewritten as a functions of H 
and H. Since q and x have a dimension of the inverse time, it is proposed to use 
the following simplest structures: 



H H 

X{t) = x^H{t) + xi^ , q{t) = qoH{t) + gi- . 



(48) 



Here Xoi Xii Qo and qi are dimensionless constants. From (21) one obtains explic- 
itly 



I{t,to 



' a{t) ' 


XO 


\ m 1 


XI 


_a{to) 




[^(^o)J 





(49) 



The function J{t, to) 

Jti 



dt 



a{t) 
a{to) 



1 XO 



m 



ixi r 



qoH{t) + qi 



H 
H 



(50) 



can be calculated explicitly for two well-known cases. These results are presented 
briefly in the following subsections 3.4.1. and 3.4.2. 



3.4.1 Power-Law Expansion 

When 



a{to) 



t \^ 
To. 



we obtain immediately the following functions: 

(7^0 - qi) 



J{t,to) 



ilXo - Xi) 



[I{t,to)-l] , I{t,to) 



f \1X0-Xl 

To) 



(51) 



(52) 



The predictions concerning the accelerated expansion depend on the sign of the 

parameter jxo — Xi- Let us consider three standard cases. 

First case: 7x0 > Xi- 

The inequahty (45) is satisfied when 



/ \ 7X0 -XI 



1^ 1 , 47(7x0 -Xi) 



2 ' ^ 4 ' D^i^qo - q^y 



(53) 



Second case: 7x0 < Xi- 
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When the discriminant is positive, i.e., I67I7X0 ~ Xi\ < D'^ilQo ~ Qifi the in- 
equahty (45) takes place, if 



1 _ 47I7X0 - Xil 
2 \ 4 L>2(^^^ _ 



f \ -|7X0-Xl| 



1^ 1 47I7X0-X1I 



2 ' ^4 L>2(^^^_^^)5 



Third case: 7x0 = Xi- 

For this special case the formulas (52) give 



I{t, to) = l, J{t, to) = (7^0 - Qi) ■ log . 



The inequahty (45) is satisfied when 

t 



log 



> 



47 



to J i.iqo-qifD'^ 



3.4.2 Exponential Expansion 

When 



a{t) 
a{to) 

the integration gives 



J(t, ^o) = - [lit, to) - 1] , I{t, to) = e>^°^°(*-*°) 
Xo 



First case: xo > 0. 

The inequahty (45) is satisfied when 

{t - to) > ^7- log 



XqHo 



1 1 I 4x0 
2^V4^2go', 



£)2 2 

Second case: Xo < 0, |xo| < 
The inequality is satisfied when 



1 


1 4x0 




4 D^qi 



< g-|xo|^^o(«-«o) ^ 



2 + ^4 



1 4|xo| 



Third case: Xo = 
The formulas (58) give 



/(t,io) = l, J{t,to)^qoHo{t-to). 
The inequality (45) is satisfied when 

4 



{t - to) > 



qiHoD^ 



(54) 



(55) 



(56) 



(57) 



(58) 



(59) 



(60) 



(61) 



(62) 
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4 Discussion 



Here the model of evolution of the self-interacting gas (fluid) Universe is pre- 
sented. The model is based on the suggestion that the expansion of the Universe 
gives rise to the specific non-equilibrium self-interaction in the kinetic system. 
The force effecting the particle can be classified as a Stokes friction force, since it 
disappears when particle co-moves with the system as a whole. The force under 
consideration can be indicated as a Langevin force, since it is proportional to 
the random scalar variable, and the particle motion can be classified as a sort 
of Brownian motion. Finally, this force can be called tidal or curvature induced 
force, since it depends on Hubble rate and its derivative, which are known to 
form the Riemann tensor, Ricci tensor, and Ricci scalar. In this sense one can 
say that we deal with curvature induced Stokes-Langevin force. 

In the subsection 3.4. the time intervals are found, during which the Universe 
expands with acceleration. The explanation of such a behaviour may be the fol- 
lowing. The self-interaction in the gas (fluid) produces the growth of the averaged 
energy and pressure in the system, in contrast to the decreasing of these param- 
eters due to the expansion. The interplay of such conflicting tendencies provides 
non-monotonic behaviour of the a function, and the model admits the existence 
of both periods of evolution: expansion with acceleration and expansion with de- 
celeration. The model under consideration requires to investigate the next step 
- the estimation of the parameters, which have been phenomenologically intro- 
duced into the force term. This task is planned to be fulfilled in the next paper 

m 
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